On exact mappings between fermionic Ising spin glass and classical spin glass models 
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We present in this paper exact analytical expressions for the thermodynamical properties and 
Green functions of a certain family of fermionic Ising spin-glass models with Hubbard interaction, 
by noticing that their Hamiltonian is a function of the number operator only. The thermodynamical 
properties are mapped to the classical Ghatak-Sherrington spin-glass model while the the Density of 
States (DoS) is related to its joint spin-field distribution. We discuss the presence of the pseudogap 
in the DoS with the help of this mapping. 
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I. INTRODUCTION 

While ferromagnetism is theoretically grounded in 
models lying between two extreme pictures, that of 
localised-spins and that of itinerant-electron theory, 
the theoretical description of spin-glass systems has 
been focused mainly on models of localised-spins, 
whose paradigms are the Edwards-Anderson (EA) and 
Sherrington-Kirkpatrick (SK) models 0, One rea- 
son for this is that experimentally most of the classic 
magnetic materials presenting spin-glass behaviour cor- 
respond to this description of localised spins, another is 
that they encapsulate many aspects of theoretical chal- 
lenge, while a third is that they relate to or emulate many 
problems of wider interest in the statistical physics of 
complex systems. 

There do exist, however, spin glass materials which are 
more appropriately described in terms of itinerant elec- 
trons, requiring models that treat magnetic and conduct- 
ing properties on the same footing. While early models 
can be found in 0, and the full problem remains to 
be tackled, our main goal here is to discuss a restricted 
class of models introduced first by Oppermann et al. @ , 
and normally referred as fermionic Ising spin-glass mod- 

ds iaaaimiiieia. 

Even within an SK-like (infinite-ranged exchange) itin- 
erant fermionic model it is a very significant challenge 
to treat conducting and magnetic properties together. 
Hence, as a first step towards their understanding, sim- 
plified models have been studied in the so-called insulat- 
ing limit [E 0, IE H E3, EH E3, O This removes the 
essentially quantum complexity of the model and allows 
a classical treatment, albeit still with interesting conse- 
quences. 

Independently of whether the insulating limit of 
fermionic Ising spin-glass models may or may not be use- 
ful to better understand real itinerant spin glasses, it is 
clear that, at least, we must fully understand the mod- 
els arising from this limit. To this end, our goal is to 
point out that not only are the fermionic Ising spin-glass 
models completely mappable to classical spin-glass mod- 
els at the level of the thermodynamics [1 1| . but also the 
densities of states (DoS), and hence the local (quantum) 
Green functions, have a classical derivation and indeed 



are given by distributions of local fields of a correspond- 
ing classical model, without the need for sophisticated 
quantum treatment. 

Hence, we can exploit all the knowledge of classical 
spin glasses to shed light on the fermionic Ising spin- 
glass models in the insulating limit. In particular, the 
existence of a pseudogap in the DoS at half-filling and 
without a Hubbard term emerges as an immediate con- 
sequence of the mapping, when account is taken of the 
well-known fact that at zero temperature the local field 
distribution of the SK spin glass has such a pseudogap. 

In turn, this implies that the observed strong correc- 
tions to the DoS due to steps in the replica symmetry 
breaking do not have a fundamentally quantum ori- 
gin. 

We also notice, en passant, the strong temperature de- 
pendence of the DoS, mirroring that of the field distribu- 
tion of the SK model. 

This paper is organised as follows: in section [H] the 
fermionic Ising spin-glass model is presented and some 
limits as a function of its parameters discussed. In section 
IIIII we map this model to the Ghatak-Sherrington model 
and express the DoS as a function of its joint spin-field 
distribution. Then in section IIVI we discuss the map- 
pings from a physical perspective and in section we 
discuss the existence of a DoS pseudogap in the light of 
the mapping. The last section of the main text is for the 
conclusions. 



II. MODEL DEFINITIONS 

Our starting point is the following model for itinerant 
electrons involving frustrated magnetic order 
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where the couplings Jij are drawn randomly and inde- 
pendently from a distribution 
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The spin and charge operators are defined by 
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(3) 



where a is and a, s are respectively the fermion creation 
and annihilation operators. The label s G {4-T} = 
{ — 1,1} indicates the spin state. The Hamiltonian 7i res t 
contains those terms that can not be expressed as a func- 
tion of the number operator only; as, for example, the 
hopping and the pair-hopping terms, viz. 



Whopping — tij a \s a jsi 

-impair \ ,pair^>f ^ ^ 
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as well as transverse spin-exchange terms. 

The class of models described by with W rC st = 
have been called fermionic Ising spin-glass models (FISG) 
|l4j and may be considered as the insulating limit of the 
larger class of itinerant models described by the Hamil- 
tonian (JIJ. Henceforth, we consider this limit. 

In the past these models have been studied using tech- 
niques of coherent fermionic states (see for instance ^3|). 
Our purpose here is to point out that quantum techniques 
are unnecessary and a classical treatment suffices. 



III. MAPPING TO CLASSICAL SPIN-GLASS 
MODELS 

For a general fermionic problem, the coherent states 
representation is a powerful technique and is likely to be 
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useful for a treatment of the full quantum Hamiltonian 
However, for the unfamiliar its use is likely to ob- 
scure a simplicity of the insulating case. Hence, here we 
proceed differently, in what we consider to be a much 
simpler way, for the Hamiltonian 
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First we note that the Hamiltonian is function of the 
number operator rii S only. We are therefore in the ideal 
position of knowing the eigenstates of the Hamiltonian 
exactly. If we define 
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then the partition function expressed in this set of states 
becomes 



.zG8)=£ e " 
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where H(n) is the Hamiltonian © with the operators 
now just numbers. This Hamiltonian is quite similar to 
such of a 3-state spin glass model. In order to make 
this similarity more apparent we express the partition 
function and the Hamiltonian as a function of the two 
new variables Si — — riii and t; = + n^. After 
doing the trace with respect to the variables r^, we end 
up with the partition function 
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Wgs(S) = - £ Wi-Bjfl?-^ log(l + e-/W), 
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with 



D = n- Tlog 1 
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and with notation S = . . . , S N ), Si £ {0, ±1}. The 
Hamiltonian (JSJ is known in classical spin glass literature 
as the Ghatak-Sherrington (GS) model [15(, a particu- 
lar case of the Blume-Emery-Griffiths-Capel spin-glass 
model |16| without biquadratic interaction. 



In the GS formulation, the new spin variables S do 
not reflect the difference between unoccupied and dou- 
bly occupied sites. It is, therefore, useful to consider the 
calculation of the expectation value of the fermion num- 
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ber operator in this formulation. By defining 

N , N 



i=l s=±l i=l 



then we can write the following relation 

_ 2(1 -p) 
n ePiu-ivC) + i + ^ ' 
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with ( • )j. the thermal average with respect a Hamilto- 
nian 7i, viz. 



= Z- 1 (/3)Tre^(---), 



(12) 



where Tr denotes the trace. Thus, half-filling (n = 1) 
corresponds to // = U/2. 

This mapping between the fermionic Ising spin-glass 
model JU and the classical SG model JHJ was already no- 
ticed and used fruitfully in |lll |. but unfortunately these 
authors appear not to have noticed that the Green func- 
tion, and the DoS which can be derived from it, can also 
be obtained simply from the classical model|2ll|. 

Indeed, instead of using the fermionic path integral 
definition for the DoS, let us start with the standard def- 
inition for the retarded Green function 

G!/(t - f) = -i6(t - if) ({a is (t),a] a ,(t')}) , (13) 



with {A, B} = AB + BA and the creation and annihi- 
lation operators given in the Heisenberg representation 
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Using the set of states |n) and after some standard ma- 
nipulations the retarded Green function takes the form 
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with 
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Fermion Ising spin-glass model without 
Hubbard interaction 



Due to the presence of | (n|ai S |m) | 2 in eq. Ijl6(l the 
only states which contribute have 



W(m) - W(n) = -sfti(S) -m, 



(17) 



with hi(S) = 'EjtyfiJijSj, the local field at site i. 
Hence, after some algebra, we can write 



dh 



u) + {fjb+ sh) /h + iS 
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Notice that all the spin dependence is in the last term. 
We proceed as as before, changing variables Si = n^—rm 
and Tj = riii + n U an( l tracing out the dependence on the 
t's. We can then re-write the above expression as 

Q'M = / dh ; ; , P?S W, t , , s (19) 



lu + (fi + sh)/h + iS 



with (h) the density of local fields at site i, viz. 



pr b w = (^[/i-/ii(s)]> WGS 



(20) 



Defining the DoS as the imaginary part of the spectral 
density function averaged over all sites, spin orientation 
and over the disorder, denoting the latter by an over line 
using the identity -^p^ — P \ — inS(x) and shifting 
the energy levels e = Hu> + /i, we obtain finally 
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with definition 



(22) 



If j7o = then we have that the distribution of fields is 
an even function, i.e. p (e) = p GS (— e), and therefore 
the expression 112 111 reveals that the DoS in the fermionic 
Ising spin-glass model without Hubbard interaction is ex- 
actly the distribution of local fields in the corresponding 
classical Ghatak-Sherrington model for any value of the 
chemical potential and temperature. 



Notice that the above expression is fully general for any 
Hamiltonian system depending on the number operator 
only. Even though we could continue with the general 
calculation quite easily, at this stage we believe it to be 
helpful first to analyse the family of fermionic Ising spin- 
glass models without Hubbard interaction, which have 
been studied extensively with coherent-state methods 9, 
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B. Fermion Ising spin-glass model with Hubbard 
interaction 

The fermionic Ising spin-glass model with Hubbard in- 
teraction was first studied in Again, due to the term 
(n| a,i S |m) | 2 , in eq. I|16(l we can replace 



H(m) - H(n) = Urns ~ sh, t (S) - fi 



(23) 
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to yield 

OfH = E / 



dh 
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where s = — s. We proceed as before and map to the GS 
model. This calculation is a bit more involved but fairly 
straightforward and after some algebra we arrive at 



#H = E E / * 

... — n i n « 



pr(r,h) 
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w + (/i + sft - 7?7) jh + iS 
(25) 



where we have introduced the joint spin-field distribution 
at site i of the GS spin-glass 



pr(T,h) = (8 8i ,r6[e-h i (S)])~ 



(26) 



with 
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From here we have that the DoS is given by 



(27) 
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with 
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In this case we have again an intimate relationship be- 
tween the DoS in the fermonic Ising spin-glass model and 
the joint spin-field distribution of the classical GS spin- 
glass model. 



IV. PHYSICAL MAPPING OF DOS TO FIELD 
DISTRIBUTIONS 

In complementation of the formal mathematical map- 
pings discussed above, in this section we describe how 
the above connection between the classical and quantum 
systems also appears naturally based solely on physical 
arguments. For the sake of simplicity we restrict the dis- 
cussion to zero temperature and Jo = 0. 

Let us consider initially that U = and fj, = 0, i.e. 
half-filling. In this case, the number of fcrmions iVf is 
equal to the number of sites N s i tcs . It is a reasonable 
ansatz, which we shall later show to be true, that in 
the ground state every site will carry a single fermion, 
whose spin can be either up or down. It is immedi- 
ately clear that this is nothing but the usual classical SK 



model. Consequently the ground state of the fermionic 
Ising spin-glass model is the same as that of the classical 
SK model and the DoS of the former system is simply 



p DoS (e)=^(-|e|). 



SK, 



where 



p SK (h) 
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(30) 



(31) 



Next let us consider further the ansatz that each site 
is singly occupied. Were a site to be unoccupied then 
clearly it would contribute no energy to the ground state. 
Neither would a doubly occupied site, since the two 
spins would both see the same effective field due to the 
other spins and they would contribute cancelling ener- 
gies. However, we also know that in the ground state 
of the SK model all local fields are finite and that spins 
are oriented to yield negative energies. Consequently it 
follows that removing a fermion from one singly occu- 
pied site and depositing it on another, already favourably 
singly occupied, site incurs two energetic penalties. Fur- 
thermore, since the distribution of local fields in the SK 
model goes to zero at zero field and any single spin cou- 
pling strength scales as ./V -1 the loss cannot be compen- 
sated by further re-adjustments on other sites. 

For \x ^ account must be taken of the fact that in the 
fermionic model some sites must be unoccupied for fi < 
or doubly occupied for fi > 0. In both cases the system 
behaves energetically as though it were a diluted classical 
SK model with spins absent on the sites of either zero or 
double occupancy in the fermionic model. Furthermore, 
the location of these holes is chosen so as to minimise 
the total ground state energy; i.e. the system behaves as 
though one has an effective Hamiltonian 

N N 

H{a, n) = - E JijWmnj ~ A*E Ui 

i<j=l i=l (-32) 

a = ±1 , n = 0, 1 , 

i.e. with two types of annealed variables, Ising spins 
(characterised by the a) and 'quasi-particles' (charac- 
terised by the rij and not to be confused with the real 
fermions of number operator s ) . We shall refer to this 
system as the anneal- diluted SK model (ADSK). The two 
chemical potentials, of cqs. and (JIJ, are related by 



I Ml 



The DoS is given by 

p DoS (£) ^ADSK He|); 



(33) 



(34) 



where p ADSK is defined analogously to p SK but with the 
sum over only the singly occupied sites and averaged over 
the ADSK Hamiltonian 



5 



It is tempting to think that the truncation of site occu- 
pation might modify the DoS of the p = case by simply 
moving the Fermi level of the p D (e) corresponding to 
the pure SK model so as to occupy only the lowest states 
up to p, but this does not take account of the loss of con- 
tribution to the fields of the unoccupied sites of (|32|) : in 
fact, computer studies of the TAP equations have shown 
that /9 D (e) now goes to zero at e = p, in a manner at 
least qualitatively similar to what happens at e = for 
the case oi p = A replica symmetric analysis of 

p close to e = p also behaves analogously to the corre- 
sponding replica symmetric study for the undiluted SK 
model near e = 0; the full replica symmetry breaking 
calculation has not yet been done explicitly. 

For p DoS (e) with |e| < \p\ it is necessary to calculate 
the local field distribution hi = J2j a ^ sites of <|32[) 

where there is no quasi-particle, so that such sites do not 
contribute to the total energy or the field or spin orien- 
tation at other sites. With this extension H34|) applies for 
all e. 

For U > the mapping of (|32[l continues to apply 
with p appropriately chosen. Again if p < U/2 there are 
(iVgite — -/Vf) unoccupied sites and for p > U/2 there are 
(Nf — A^ite) doubly occupied fermion sites; consequently 
in both cases there are \ (N{ — A s i tc )| sites without quasi- 
particles. p is given by |22j| 



r /i, fi < u/2 

' \U - p, p>U/2 



(35) 



and p Dos (e) is given by 



p DoS (e) = { 



9 ADSK (-|e|) 



_,ADSK 



e < 

< e < U (36) 



(-\e-U\) e>0. 



These mappings may be related to those of the pre- 
vious section by noting that Hamiltonian (|32|) is also 
another way of writing the GS model of JSJ with p = 
D-T\n2 



V. EXISTENCE OF A PSEUDOGAP IN 
FERMIONIC ISING SPIN GLASSES 

Having demonstrated the mapping between DoS and 
distribution of fields, we can draw some conclusions and 
speculate about the nature of the pseudogap in the DoS 
at the Fermi energy [l^. First of all, let us notice that 
the nature of the pseudogap and of the strong corrections 
of different steps of replica symmetry breaking, cannot be 
but of classical origin. The strong corrections to the DoS 
from RSB corrections found in fermionic Ising spin-glass 
models are common in classical spin glasses. 

In particular is has been shown already for some 
time from oo— RSB calculations 0, and Thouless- 
Anderson-Palmer equations [l^, that the field distribu- 
tion at zero temperature in the SK model vanishes with 



h -> as p SK {h) = a\h\ (See also (U). This therefore 
predicts the existence of a pseudogap for the fermionic 
Ising-spin glass model without Hubbard term and at half- 
filling 111. 

The distribution of fields for fermonic SK models with- 
out Hubbard interaction for p ^ was studied numeri- 
cally using a TAP approach in [l(| for p > 0, indicating 
that it seems to vanish, again quasi-linearly, at h = ±/i, 
implying that the DoS also presents pseudogaps. 

We might note also that the field distribution is very 
temperature-dependent and the pseudogap becomes 
filled in as the temperature rises; so therefore (and 
unusually) the density of fermionic states will mirror 
this strong temperature-dependence. 



VI. CONCLUSIONS 

In this paper we have shown that the fermionic Ising 
spin-glass model (with SK-like interactions) is mappable 
to the classical GS spin-glass model not only at the level 
of the free energy but also the DoS, the latter being 
given exactly by the local field distribution of the GS 
model in the case without Hubbard interaction and ob- 
tainable from it when a Hubbard term is present. By 
using known results from spin glass models we can show 
the existence of a pseudogap, from full RSB and TAP 
approachs. It should be noted that the pseudogap and 
strong corrections in the different steps of RSB are purely 
classical effects and not due to quantum fluctuations. It 
would be interesting to see how this picture changes when 
the Hamiltonian H rcs t is switched on and also when one 
passes to a more realistic (but also more difficult to solve 
and currently still controversial) model with short-range 
interactions. 

Above we have concentrated on single-fermion Green 
functions and their averages. A similar procedure to that 
outlined in Section III can be applied to higher-order 
Green functions, mapping into higher order field distri- 
butions, and for averages of products of Green functions. 

Finally, we note that the GS model has a first-order 
phase transition at a critical negative- valued D C (T) and 
beneath a tri-critical temperature T3, between a mag- 
netic state (| {Si) I ^ 0) and a non-magnetic (| (Si) | = 0) 
solution. This reflects in the FISG to a critical p c for 
magnetic breakdown [23j . 
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